Detection of driver mutations and genomic signatures in endometrial cancers using artificial intelligence algorithms.

Analyzed endometrial cancer (EC) genomes have allowed for the identification of molecular signatures, which enable the classification, and sometimes prognostication, of these cancers. Artificial intelligence algorithms have facilitated the partitioning of mutations into driver and passenger based on a variety of parameters, including gene function and frequency of mutation. Here, we undertook an evaluation of EC cancer genomes deposited on the Catalogue of Somatic Mutations in Cancers (COSMIC), with the goal to classify all mutations as either driver or passenger. Our analysis showed that approximately 2.5% of all mutations are driver and cause cellular transformation and immortalization. We also characterized nucleotide level mutation signatures, gross chromosomal re-arrangements, and gene expression profiles. We observed that endometrial cancers show distinct nucleotide substitution and chromosomal re-arrangement signatures compared to other cancers. We also identified high expression levels of the CLDN18 claudin gene, which is involved in growth, survival, metastasis and proliferation. We then used in silico protein structure analysis to examine the effect of certain previously uncharacterized driver mutations on protein structure. We found that certain mutations in CTNNB1 and TP53 increase protein stability, which may contribute to cellular transformation. While our analysis retrieved previously classified mutations and genomic alterations, which is to be expected, this study also identified new signatures. Additionally, we show that artificial intelligence algorithms can be effectively leveraged to accurately predict key drivers of cancer. This analysis will expand our understanding of ECs and improve the molecular toolbox for classification, diagnosis, or potential treatment of these cancers.

F-box DNA Helicase 1 (FBH1) Contributes to the Destabilization of DNA Damage Repair Machinery in Human Cancers.

Homologous recombination (HR) is the major mechanism of rescue of stalled replication forks or repair of DNA double-strand breaks (DSBs) during S phase or mitosis. In human cells, HR is facilitated by the BRCA2-BRCA1-PALB2 module, which loads the RAD51 recombinase onto a resected single-stranded DNA end to initiate repair. Although the process is essential for error-free repair, unrestrained HR can cause chromosomal rearrangements and genome instability. F-box DNA Helicase 1 (FBH1) antagonizes the role of BRCA2-BRCA1-PALB2 to restrict hyper-recombination and prevent genome instability. Here, we analyzed reported FBH1 mutations in cancer cells using the Catalogue of Somatic Mutations in Cancers (COSMIC) to understand how they interact with the BRCA2-BRCA1-PALB2. Consistent with previous results from yeast, we find that FBH1 mutations co-occur with BRCA2 mutations and to some degree BRCA1 and PALB2. We also describe some co-occurring mutations with RAD52, the accessory RAD51 loader and facilitator of single-strand annealing, which is independent of RAD51. In silico modeling was used to investigate the role of key FBH1 mutations on protein function, and a Q650K mutation was found to destabilize the protein structure. Taken together, this work highlights how mutations in several DNA damage repair genes contribute to cellular transformation and immortalization.

Protein Arginine Methyltransferase 5 (PRMT5) Mutations in Cancer Cells.

Arginine methylation is a form of posttranslational modification that regulates many cellular functions such as development, DNA damage repair, inflammatory response, splicing, and signal transduction, among others. Protein arginine methyltransferase 5 (PRMT5) is one of nine identified methyltransferases, and it can methylate both histone and non-histone targets. It has pleiotropic functions, including recruitment of repair machinery to a chromosomal DNA double strand break (DSB) and coordinating the interplay between repair and checkpoint activation. Thus, PRMT5 has been actively studied as a cancer treatment target, and small molecule inhibitors of its enzymatic activity have already been developed. In this report, we analyzed all reported PRMT5 mutations appearing in cancer cells using data from the Catalogue of Somatic Mutations in Cancers (COSMIC). Our goal is to classify mutations as either drivers or passengers to understand which ones are likely to promote cellular transformation. Using gold standard artificial intelligence algorithms, we uncovered several key driver mutations in the active site of the enzyme (D306H, L315P, and N318K). In silico protein modeling shows that these mutations may affect the affinity of PRMT5 for S-adenosylmethionine (SAM), which is required as a methyl donor. Electrostatic analysis of the enzyme active site shows that one of these mutations creates a tunnel in the vicinity of the SAM binding site, which may allow interfering molecules to enter the enzyme active site and decrease its activity. We also identified several non-coding mutations that appear to affect PRMT5 splicing. Our analyses provide insights into the role of PRMT5 mutations in cancer cells. Additionally, since PRMT5 single molecule inhibitors have already been developed, this work may uncover future directions in how mutations can affect targeted inhibition.

In silico analysis of a SLC6A4 G100V mutation in lung cancers.

SLC6A4 is a serotonin re-uptake transporter which has been a target for anti-depressant therapies but recently some mutations have been described in cancer cells. Here, we characterize mutations in SLC6A4 that appear in cancer cells. We employed several validated computational and artificial intelligence algorithms to characterize the mutations. We identified a previously uncharacterized G100V mutation in lung cancers. In sillico structural analysis reveals that this mutation may affect SLC6A4 ligand binding and subsequently its function. We also identified several other mutations that may affect the structure of the protein. This preliminary analysis highlights the role of SLC6A4 in human cancers.

Pan-cancer analysis of co-occurring mutations in RAD52 and the BRCA1-BRCA2-PALB2 axis in human cancers.

In human cells homologous recombination (HR) is critical for repair of DNA double strand breaks (DSBs) and rescue of stalled or collapsed replication forks. HR is facilitated by RAD51 which is loaded onto DNA by either BRCA2-BRCA1-PALB2 or RAD52. In human culture cells, double-knockdowns of RAD52 and genes in the BRCA1-BRCA2-PALB2 axis are lethal. Mutations in BRCA2, BRCA1 or PALB2 significantly impairs error free HR as RAD51 loading relies on RAD52 which is not as proficient as BRCA2-BRCA1-PALB2. RAD52 also facilitates Single Strand Annealing (SSA) that produces intra-chromosomal deletions. Some RAD52 mutations that affect the SSA function or decrease RAD52 association with DNA can suppress certain BRCA2 associated phenotypes in breast cancers. In this report we did a pan-cancer analysis using data reported on the Catalogue of Somatic Mutations in Cancers (COSMIC) to identify double mutants between RAD52 and BRCA1, BRCA2 or PALB2 that occur in cancer cells. We find that co-occurring mutations are likely in certain cancer tissues but not others. However, all mutations occur in a heterozygous state. Further, using computational and machine learning tools we identified only a handful of pathogenic or driver mutations predicted to significantly affect the function of the proteins. This supports previous findings that co-inactivation of RAD52 with any members of the BRCA2-BRCA1-PALB2 axis is lethal. Molecular modeling also revealed that pathogenic RAD52 mutations co-occurring with mutations in BRCA2-BRCA1-PALB2 axis are either expected to attenuate its SSA function or its interaction with DNA. This study extends previous breast cancer findings to other cancer types and shows that co-occurring mutations likely destabilize HR by similar mechanisms as in breast cancers.

Distribution of copy number variations and rearrangement endpoints in human cancers with a review of literature.

Copy number variations (CNVs) which include deletions, duplications, inversions, translocations, and other forms of chromosomal re-arrangements are common to human cancers. In this report we investigated the pattern of these variations with the goal of understanding whether there exist specific cancer signatures. We used re-arrangement endpoint data deposited on the Catalogue of Somatic Mutations in Cancers (COSMIC) for our analysis. Indeed, we find that human cancers are characterized by specific patterns of chromosome rearrangements endpoints which in turn result in cancer specific CNVs. A review of the literature reveals tissue specific mutations which either drive these CNVs or appear as a consequence of CNVs because they confer an advantage to the cancer cell. We also identify several rearrangement endpoints hotspots that were not previously reported. Our analysis suggests that in addition to local chromosomal architecture, CNVs are driven by the internal cellular or nuclear physiology of each cancer tissue.

Arginine Depletion in Human Cancers.

Arginine is encoded by six different codons. Base pair changes in any of these codons can have a broad spectrum of effects including substitutions to twelve different amino acids, eighteen synonymous changes, and two stop codons. Four amino acids (histidine, cysteine, glutamine, and tryptophan) account for over 75% of amino acid substitutions of arginine. This suggests that a mutational bias, or "purifying selection", mechanism is at work. This bias appears to be driven by C > T and G > A transitions in four of the six arginine codons, a signature that is universal and independent of cancer tissue of origin or histology. Here, we provide a review of the available literature and reanalyze publicly available data from the Catalogue of Somatic Mutations in Cancer (COSMIC). Our analysis identifies several genes with an arginine substitution bias. These include known factors such as IDH1, as well as previously unreported genes, including four cancer driver genes (FGFR3, PPP6C, MAX, GNAQ). We propose that base pair substitution bias and amino acid physiology both play a role in purifying selection. This model may explain the documented arginine substitution bias in cancers.

A Genome-Wide Profiling of Glioma Patients with an IDH1 Mutation Using the Catalogue of Somatic Mutations in Cancer Database.

Gliomas are differentiated into two major disease subtypes, astrocytoma or oligodendroglioma, which are then characterized as either IDH (isocitrate dehydrogenase)-wild type or IDH-mutant due to the dramatic differences in prognosis and overall survival. Here, we investigated the genetic background of IDH1-mutant gliomas using the Catalogue of Somatic Mutations in Cancer (COSMIC) database. In astrocytoma patients, we found that IDH1 is often co-mutated with TP53, ATRX, AMBRA1, PREX1, and NOTCH1, but not CHEK2, EGFR, PTEN, or the zinc finger transcription factor ZNF429. The majority of the mutations observed in these genes were further confirmed to be either drivers or pathogenic by the Cancer-Related Analysis of Variants Toolkit (CRAVAT). Gene expression analysis showed down-regulation of DRG2 and MSN expression, both of which promote cell proliferation and invasion. There was also significant over-expression of genes such as NDRG3 and KCNB1 in IDH1-mutant astrocytoma patients. We conclude that IDH1-mutant glioma is characterized by significant genetic changes that could contribute to a better prognosis in glioma patients.

Comprehensive Genetic Analysis of DGAT2 Mutations and Gene Expression Patterns in Human Cancers.

DGAT2 is a transmembrane protein encoded by the DGAT2 gene that functions in lipid metabolism, triacylglycerol synthesis, and lipid droplet regulation. Cancer cells exhibit altered lipid metabolism and mutations in DGAT2 may contribute to this state. Using data from the Catalogue of Somatic Mutations in Cancer (COSMIC), we analyzed all cancer genetic DGAT2 alterations, including mutations, copy number variations and gene expression. We find that several DGAT2 mutations fall within the catalytic site of the enzyme. Using the Variant Effect Scoring Tool (VEST), we identify multiple mutations with a high likelihood of contributing to cellular transformation. We also found that D222V is a mutation hotspot neighboring a previously discovered Y223H mutation that causes Axonal Charcot-Marie-Tooth disease. Remarkably, Y223H has not been detected in cancers, suggesting that it is inhibitory to cancer progression. We also identify several single nucleotide polymorphisms (SNP) with high VEST scores, indicating that certain alleles in human populations have a pathogenic predisposition. Most mutations do not correlate with a change in gene expression, nor is gene expression dependent on high allele copy number. However, we did identify eight alleles with high expression levels, suggesting that at least in certain cases, the excess DGAT2 gene product is not inhibitory to cellular proliferation. This work uncovers unknown functions of DGAT2 in cancers and suggests that its role may be more complex than previously appreciated.

Schizosaccharomyces pombe KAT5 contributes to resection and repair of a DNA double-strand break.

Chromatin remodeling is essential for effective repair of a DNA double-strand break (DSB). KAT5 (Schizosaccharomyces pombe Mst1, human TIP60) is a MYST family histone acetyltransferase conserved from yeast to humans that coordinates various DNA damage response activities at a DNA DSB, including histone remodeling and activation of the DNA damage checkpoint. In S. pombe, mutations in mst1+ causes sensitivity to DNA damaging drugs. Here we show that Mst1 is recruited to DSBs. Mutation of mst1+ disrupts recruitment of repair proteins and delays resection. These defects are partially rescued by deletion of pku70, which has been previously shown to antagonize repair by homologous recombination (HR). These phenotypes of mst1 are similar to pht1-4KR, a nonacetylatable form of histone variant H2A.Z, which has been proposed to affect resection. Our data suggest that Mst1 functions to direct repair of DSBs toward HR pathways by modulating resection at the DSB.

Mutation Spectra of the MRN (MRE11, RAD50, NBS1/NBN) Break Sensor in Cancer Cells.

The MRN complex (MRE11, RAD50, NBS1/NBN) is a DNA double strand break sensor in eukaryotes. The complex directly participates in, or coordinates, several activities at the break such as DNA resection, activation of the DNA damage checkpoint, chromatin remodeling and recruitment of the repair machinery. Mutations in components of the MRN complex have been described in cancer cells for several decades. Using the Catalogue of Somatic Mutations in Cancer (COSMIC) database, we characterized all the reported MRN mutations. This analysis revealed several hotspot frameshift mutations in all three genes that introduce premature stop codons and truncate large regions of the C-termini. We also found through evolutionary analyses that COSMIC mutations are enriched in conserved residues of NBS1/NBN and RAD50 but not in MRE11. Given that all three genes are important to carcinogenesis, we propose these differential enrichment patterns may reflect a more severe pleiotropic role for MRE11.

Identifying Cancer-Relevant Mutations in the DLC START Domain Using Evolutionary and Structure-Function Analyses.

Rho GTPase signaling promotes proliferation, invasion, and metastasis in a broad spectrum of cancers. Rho GTPase activity is regulated by the deleted in liver cancer (DLC) family of bona fide tumor suppressors which directly inactivate Rho GTPases by stimulating GTP hydrolysis. In addition to a RhoGAP domain, DLC proteins contain a StAR-related lipid transfer (START) domain. START domains in other organisms bind hydrophobic small molecules and can regulate interacting partners or co-occurring domains through a variety of mechanisms. In the case of DLC proteins, their START domain appears to contribute to tumor suppressive activity. However, the nature of this START-directed mechanism, as well as the identities of relevant functional residues, remain virtually unknown. Using the Catalogue of Somatic Mutations in Cancer (COSMIC) dataset and evolutionary and structure-function analyses, we identify several conserved residues likely to be required for START-directed regulation of DLC-1 and DLC-2 tumor-suppressive capabilities. This pan-cancer analysis shows that conserved residues of both START domains are highly overrepresented in cancer cells from a wide range tissues. Interestingly, in DLC-1 and DLC-2, three of these residues form multiple interactions at the tertiary structural level. Furthermore, mutation of any of these residues is predicted to disrupt interactions and thus destabilize the START domain. As such, these mutations would not have emerged from traditional hotspot scans of COSMIC. We propose that evolutionary and structure-function analyses are an underutilized strategy which could be used to unmask cancer-relevant mutations within COSMIC. Our data also suggest DLC-1 and DLC-2 as high-priority candidates for development of novel therapeutics that target their START domain.

Comprehensive Analysis of MEN1 Mutations and Their Role in Cancer.

MENIN is a scaffold protein encoded by the MEN1 gene that functions in multiple biological processes, including cell proliferation, migration, gene expression, and DNA damage repair. MEN1 is a tumor suppressor gene, and mutations that disrupts MEN1 function are common to many tumor types. Mutations within MEN1 may also be inherited (germline). Many of these inherited mutations are associated with a number of pathogenic syndromes of the parathyroid and pancreas, and some also predispose patients to hyperplasia. In this study, we cataloged the reported germline mutations from the ClinVar database and compared them with the somatic mutations detected in cancers from the Catalogue of Somatic Mutations in Cancer (COSMIC) database. We then used statistical software to determine the probability of mutations being pathogenic or driver. Our data show that many confirmed germline mutations do not appear in tumor samples. Thus, most mutations that disable MEN1 function in tumors are somatic in nature. Furthermore, of the germline mutations that do appear in tumors, only a fraction has the potential to be pathogenic or driver mutations.

Secondary Resistant Mutations to Small Molecule Inhibitors in Cancer Cells.

Secondary resistant mutations in cancer cells arise in response to certain small molecule inhibitors. These mutations inevitably cause recurrence and often progression to a more aggressive form. Resistant mutations may manifest in various forms. For example, some mutations decrease or abrogate the affinity of the drug for the protein. Others restore the function of the enzyme even in the presence of the inhibitor. In some cases, resistance is acquired through activation of a parallel pathway which bypasses the function of the drug targeted pathway. The Catalogue of Somatic Mutations in Cancer (COSMIC) produced a compendium of resistant mutations to small molecule inhibitors reported in the literature. Here, we build on these data and provide a comprehensive review of resistant mutations in cancers. We also discuss mechanistic parallels of resistance.

Schizosaccharomyces pombe Assays to Study Mitotic Recombination Outcomes.

The fission yeast-Schizosaccharomyces pombe-has emerged as a powerful tractable system for studying DNA damage repair. Over the last few decades, several powerful in vivo genetic assays have been developed to study outcomes of mitotic recombination, the major repair mechanism of DNA double strand breaks and stalled or collapsed DNA replication forks. These assays have significantly increased our understanding of the molecular mechanisms underlying the DNA damage response pathways. Here, we review the assays that have been developed in fission yeast to study mitotic recombination.

An Assay to Study Intra-Chromosomal Deletions in Yeast.

An accurate DNA damage response pathway is critical for the repair of DNA double-strand breaks. Repair may occur by homologous recombination, of which many different sub-pathways have been identified. Some recombination pathways are conservative, meaning that the chromosome sequences are preserved, and others are non-conservative, leading to some alteration of the DNA sequence. We describe an in vivo genetic assay to study non-conservative intra-chromosomal deletions at regions of non-tandem direct repeats in Schizosaccharomyces pombe. This assay can be used to study both spontaneous breaks arising during DNA replication and induced double-strand breaks created with the S. cerevisiae homothallic endonuclease (HO). The preliminary genetic validation of this assay shows that spontaneous breaks require rad52+ but not rad51+, while induced breaks require both genes, in agreement with previous studies. This assay will be useful in the field of DNA damage repair for studying mechanisms of intra-chromosomal deletions.

Development of the SapI/AarI Incision Mediated Plasmid Editing Method.

Plasmid engineering and molecular cloning is a virtually ubiquitous tool in biology. Although various methods have been developed for ligating DNA molecules or targeted mutagenesis of plasmids, each has its limitations. Many of the commonly used laboratory strategies are inefficient, while commercially available kits are quite costly and often specialized for highly specific circumstances. Here, we describe the SapI/AarI incision mediated plasmid editing (SIMPLE) method, which allows users to perform site-directed mutagenesis, deletions, and even short insertions into any plasmid in a single PCR reaction, using just one restriction enzyme. In addition, the SIMPLE method can be adapted to insert any sized DNA fragment into a vector using a two-step PCR approach, and can be used to ligate any number of DNA fragments with non-compatible ends in the specific order desired. The SIMPLE method provides researches an efficient and powerful tool with a broad range of applications for molecular cloning.

Replication fork stability is essential for the maintenance of centromere integrity in the absence of heterochromatin.

The centromere of many eukaryotes contains highly repetitive sequences marked by methylation of histone H3K9 by Clr4(KMT1). This recruits multiple heterochromatin proteins, including Swi6 and Chp1, to form a rigid centromere and ensure accurate chromosome segregation. In the absence of heterochromatin, cells show an increased rate of recombination in the centromere, as well as chromosome loss. These defects are severely aggravated by loss of replication fork stability. Thus, heterochromatin proteins and replication fork protection mechanisms work in concert to prevent abnormal recombination, preserve centromere integrity, and ensure faithful chromosome segregation.

TEN1 is essential for CDC13-mediated telomere capping.

Telomere binding proteins protect chromosome ends from degradation and mask chromosome termini from checkpoint surveillance. In Saccharomyces cerevisiae, Cdc13 binds single-stranded G-rich telomere repeats, maintaining telomere integrity and length. Two additional proteins, Ten1 and Stn1, interact with Cdc13 but their contributions to telomere integrity are not well defined. Ten1 is known to prevent accumulation of aberrant single-stranded telomere DNA; whether this results from defective end protection or defective telomere replication is unclear. Here we report our analysis of a new group of ten1 temperature-sensitive (ts) mutants. At permissive temperatures, ten1-ts strains display greatly elongated telomeres. After shift to nonpermissive conditions, however, ten1-ts mutants accumulate extensive telomeric single-stranded DNA. Cdk1 activity is required to generate these single-stranded regions, and deleting the EXO1 nuclease partially suppresses ten1-ts growth defects. This is similar to cdc13-1 mutants, suggesting ten1-ts strains are defective for end protection. Moreover, like Cdc13, our analysis reveals Ten1 promotes de novo telomere addition. Interestingly, in ten1-ts strains at high temperatures, telomeric single-stranded DNA and Rad52-YFP repair foci are strongly induced despite Cdc13 remaining associated with telomeres, revealing Cdc13 telomere binding is not sufficient for end protection. Finally, unlike cdc13-1 mutants, ten1-ts strains display strong synthetic interactions with mutations in the POLalpha complex. These results emphasize that Cdc13 relies on Ten1 to execute its essential function, but leave open the possibility that Ten1 has a Cdc13-independent role in DNA replication.